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We perform a generalized analysis of data from WIMP search experiments for 
\ point-like WIMPs of arbitrary spin and general Lorenz-invariant WIMP-nucleus in- 

teraction. We show that in the non-relativistic limit only spin-independent (SI) 
and spin-dependent (SD) WIMP-nucleon interactions survive, which can be param- 
eterized by only five independent parameters. We explore this five-dimensional pa- 
rameter space to determine whether the annual modulation observed in the DAMA 

> ■ 

■ experiment can be consistent with all other experiments. The pure SI interaction 
OO ■ 

, is ruled out except for very small region of parameter space with the WIMP mass 

I close to 50 GeV and the ratio of the WIMP-neutron to WIMP-proton SI couplings 

■ —0.77 < fn/ fp ^ —0.75. For the predominantly SD interaction, we find an upper 
Ph! limit to the WIMP mass of about 18 GeV, which can only be weakened if the con- 
Q \ straint stemming from null searches for energetic neutrinos from WIMP annihilation 



the Sun is evaded. None of the regions of the parameter space that can reconcile all 



^ . WIMP search results can be easily accommodated in the minimal supersymmetric 

H ; 

5^ , extension of the standard model. 

PACS numbers: 95.35.+d, 



I. INTRODUCTION 

A weakly-interacting massive particle (WIMP) is perhaps the most natural candidate 
for the dark matter that permeates our Galactic halo . Simply stated, if new physics 
at the electroweak scale introduces a neutral stable particle, then that particle has a cos- 
mo logical density comparable to that contributed by halo dark matter. The most widely 
studied possibility for the WIMP is the neutralino, a linear combination of the supersym- 
metric partners of the photon and and Higgs boson, the lightest superpartner (LSP) in 
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many minimal supersymmetric extensions of the standard model (MSSMs). However, other 
possibilities include the sneutrino (the neutrino's superpartner), heavy neutrinos, particles 
in non-minimal supersymmetric models, and Kaluza-Klein modes in models with universal 
extra dimensions; there may be others as well. 

If such WIMPs exist in the Galactic halo, they may be detected directly via detection of 
the 0(30 keV) recoil energy imparted to a nucleus in a low-background detector when a halo 
WIMP elastically collides with the nucleus j3|,|^. They may also be detected indirectly via 
observation of energetic neutrinos (i.e., E,^ ^ 10 GeV) produced by annihilation of WIMPs 
that have accumulated in the Sun and/or Earth [5]. A third possibility is observation of 
cosmic-ray antiprotons, positrons, or gamma rays produced by annihilation of dark-matter 
particles in the Galactic halo, although the certainty with which we can predict the fluxes of 
such cosmic rays is limited by considerable uncertainties in the distribution of dark matter 
in the halo as well as in cosmic-ray propagation. 

Of the several direct dark-matter-detection experiments underway, one — the DAMA 
collaboration — has reported an annual modulation in their Nal detector that they attribute 
to the difference in the WIMP flux incident on the Earth between the summer and win- 
ter 6j. Since their announcement, several other experiments have reported 
null results that they argue conflict with the DAMA detection. In order to come to this 
conclusion, however, it is assumed that the WIMP has only scalar (spin-independent; SI) 
interactions with nucleons, and with equal interaction strength with neutrons and protons. 
Although such interactions are favored in most of the supersymmetric models that appear 
in the literature, they are by no means universal to the MSSM, nor to WIMP models more 
generally. Ref. jll| studied the possibility that the DAMA results might be due to a WIMP 
with an axial- vector (spin-dependent; SD) coupling to either protons or to neutrons, and 
showed that the former possibility conflicts with null searches for energetic neutrinos from 
the Sun, while the latter is inconsistent with other direct-detection experiments. However, 
although the realm of possibilities that have been explored has been expanded, it is still not 
universal. 

In this paper we expand further on prior analyses of WIMP-detection experiments and 
explore whether the DAMA modulation may be consistent with null searches in a more 
generalized parameter space. We suppose that the WIMP is a pointlike particle, but with 
arbitrary spin. We consider particles that are either self-charge-conjugate, or that have 
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an antiparticle, either with or without a particle-antiparticle asymmetry. We begin in the 
next Section by showing that from the point of view of dark-matter detection, the Dirac 
structure of the interaction Lagrangian of all such particles is such that the interactions 
relevant for dark-matter detection can be parameterized by only five quantities: the particle 
mass, SI couplings to neutrons and protons, and SD couplings to protons and neutrons. In 
Section III we review the cross sections. Section IV discusses the experimental results we 
use from the DAMA Nal experiment, the DAMA Xe experiment, EDELWEISS, ZEPLIN-I, 
and energetic-neutrino searches Q].^ In Section V we search the five-dimensional parameter 
space for regions that are consistent with all of the experimental results. The results are 
summarized in Section VI, the Conclusions. The bottom line is that there are regions in 
this five-dimensional parameter space for both SD and SI interactions that are consistent 
with all current experimental results, but these regions are unlikely to occur in any MSSM 
and do not apply to any of the other WIMP candidates that have appeared in the recent 
literature. 



II. GENERALIZED WIMP-NUCLEON INTERACTION 



WIMP searches have usually been interpreted within the framework of the minimal su- 
persymmetric standard model (MSSM) because it contains an excellent WIMP candidate, 
the lightest neutralino. Interactions of the neutralino with quarks (and ultim ately with 
nucleons and nuclei) have been extensively studied in the literature (see e.g. I3)- For 



slow-moving Majorana neutralinos, the neutralino-nucleon interaction has only two terms. 



<Six) = fppp + fnun , 



(1) 



where .s the Ferm. constant, and and are, respectively, the proton and neutron 
coupling constants given in Ref. The term containing S{x) gives rise to the spin- 



^ We do not explicitly consider the constraints from the CDMS experiment 7] as they are surpassed by the 
constraints from the EDELWEISS experiment Q in the region of interest; both CDMS and EDELWEISS 
used cryogenic Ge detectors. 

^ We absorbed the coupling constants into the definitions of and S to make this simple notation consistent 
with the convention adopted in Ref. 0. 
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independent (SI) interaction, whereas the J^{x) term is responsible for the spin- dependent 
(SD) interaction. 

Although Eq. (0) was derived assuming that WIMPs are slow-moving Majorana spin-1/2 
fermions within the MSSM, it can be generalized to particles of arbitrary spin, indepen- 
dent of whether they are self-charge-conjugate or not. Consider a general Lorentz-invariant 
interaction for the WIMP-nucleon coupling. 



-iV75iV 



Here, S-^ and V-^ are linear combinations of scalar and pseudoscalar operators built from 
the X field, and are the corresponding vector and axial-vector operators, and T^^ 
and are the tensor and pseudotensor operators. In the Weyl representation, the Dirac 
structures appearing in Eq. ^ are (Latin indices indicate spatial components), 

/ 



+ [T^^ + I)f ) {GtNa,,N + GaNa^^^^N) . (2) 




7° = h 7^ = I h 75 




-1 
1 



a^^ = 1{7V-7V}, (3) 

where a* are Pauli matrices. In order to use this interaction to calculate matrix elements 
for WIMP-nucleus scattering, one must sum over all nucleons in the target nucleus. 

The WIMP-nucleon interaction is determined by the underlying WIMP-quark interaction, 
and to obtain Eq. Q from the underlying interaction one must take its matrix elements 
between one-nucleon states. In the limit |g| <^ Mn relevant for dark- matter detection, one 
can do this by simply replacing the quark fields with the nucleon fields and rescaling the 
corresponding coupling constants. For example, gy{N\q'~^fj_q\N) = GyUN'JfiUN, where Un is 
the nucleon spinor. This procedure holds for all operators except for the quark axial-vector 
current whose matrix element has a pole: 



{N\J^,\N) = N (gAYl' + gp{q')q^l,) N 



2n 2Mjv , . 

gpii ) = gA— — 2 ' (^) 

— q^ 

where the second line implements the PCAC hypothesis, and m,r is the pion mass. This 

argument justifies the appearance of the particular structure multiplying Ga in Eq. The 
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second term in the first line of Eq. ^ becomes important if |g| > [l3|]. For simplicity, we 
do not explicitly keep the pion pole term in the following. However, it should be understood 
that this term always appears in the form given by the last term in the second line of Eq. 

Except for the pion pole term, Eq. does not contain operators with explicit factors of 
momentum transfer q. Such terms are considered subleading. As will be shown below, such 
operators are suppressed by powers of q/Mp or q/M^. It is easy to show that in the frame 
where the nucleus is initially at rest the components of q are related by 2goM„„c + Q'^ = 0, 
which also shows that go ^ I*?!- Therefore, we find. 



where A < 100 is the atomic number of the nucleus and go is the recoil energy, which we 
assume not to significantly exceed 100 keV. 

Let us look at each of the one-nucleon operators appearing in Eq. and determine 
which of them have non- vanishing nuclear matrix elements in the non-relativistic limit. 

Scalar operators. The scalar operator NN simply counts nucleons; it obviously survives 
in the non-relativistic limit. On the other hand, nuclear matrix elements of the pseudoscalar 
operator Vn = J2k ^k{x + Vkj-^^Nk^x + r^) are suppressed: {f\VN\i) ~ \'i\lMp. Here, x and 
Tk correspond, respectively, to the center of mass of the nucleus and to the position of the 
fcth nucleon relative to the center of mass. Since our discussion is not affected by the fact 
that the nucleus has a finite size we can consider it to be a point particle and set all 
to zero. In a more accurate treatment one can account for finite size R by introducing a 
form-factor F{\q\R). Defining an "axial-vector current" operator A'i^{x) = VN{x)dx^ as 
an intermediate step we obtain: 



only axial vector describing a state of the nucleus that does not vanish as g ^ 0. The 
scalar constant P has the dimension of inverse mass. The largest quantity that has this 
property and is not singular in the limit g ^ is the inverse mass of the nucleon. Therefore, 
P ~ 1/Mp, and matrix elements of Vn are suppressed according to Eq. (0). Another way to 
see that P ~ 1/Mp (and, e.g., not 1/Mnuc) is to first consider the one-nucleon operator N'j^N 





^ ^g,e^«--'P(/|J^|z) + 0(g2) , 



(6) 



where Jjv is the nuclear spin operator. In the last step we used the fact that spin is the 
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using free Dirac spinors for the nucleon fields. Using free spinors is a good approximation 
because the binding energy of a nucleon in a nucleus is much smaller than the nucleon mass, 
which means that nucleons are only slightly off-shell. A direct calculation yields, 

N{p + q)^,N{p) = ^ya,r^ , (7) 

where 77, rj' are Weyl spinors characterizing the states of the nucleon with the initial momen- 
tum p, and (Ji is the vector of Pauli matrices. In a nucleus one obtains 

i^(/lE^IWK)~i^(/I^^K), (8) 

IVlp ^ IVlp 

in agreement with Eq. (jHl). 

Vector operators. For Vj^ = Xlfc Nk'-f^N^ we must have: 

(/IV^K) = V+{pf+p,r + V^ipf-pd'- , (9) 

where and Vj^ are constants. Here, Pi = {Mnuc, 0) and Pf = q + Pi- Just like for Vn, we 
may conclude on dimensional grounds that V^,V^ ~ ^/Mp. Then the time component of 
Vn is of order unity, whereas the spatial components are suppressed as ~ <t l^p- The 
time component multiplies V^^^o and Ay^^ in Eq. which transform, respectively, as 
a scalar and a pseudoscalar under the extended rotational group. Therefore, they can be 
effectively absorbed into and V-^ in the non-relativistic case. 

For the axial-vector operator ^^^5 = J2k ^kl^l5^k the situation is reversed. The time 
component transforms as a pseudoscalar under the extended rotational group and its ma- 
trix elements are suppressed for the same reason as those of Vn- The spatial components 
transform like a pseudovector leading to 

{f\A'^,\z) = AN{fm^) , (10) 

where An is a dimensionless constant, which can in general be of order unity. 

Tensor operators. Let us first consider the operator TJ^*^ = J2k ^kO''^'^ Nk- Nonzero 
components of a^'^ are cr°* and i,j = 1,2,3. Under the extended rotational group, a^* 
transforms as a polar vector. Therefore, in analogy to Eq. (j^I), 

= niPf + p^y + Tj^iPf - p^r , (n) 

which is suppressed similarly to V]^ in Eq. Q. To analyze TJ^ consider its dual = 
e^^^T^jk = Jlk^k^^^^cTjkNk = '2J2k ^kC^Nk, where e*-'^ is the three-dimensional Levi-Civita 
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tensor. In the last step, we used the property ajk = e^^^ci, which can be derived from Eq. Q. 
We obtain 

{f\T^\i) = U^'r^, = e'^'{f\Y.Nka'Nk\z) = T^e'^'{f\JNA^) , (12) 
^ k 

where T/v is a dimensionless constant. Therefore, survives in the non-relativistic hmit. 

The operator V'^ = Y^k^kO'^^l^N^ can be expressed as {i/2)e^^^P^Tjqpx, which follows 
from the identity = {i/2)e^''^^apx. Since only survives in the non-relativistic limit 

we conclude that are the only non- vanishing components of in this limit: 

v'j^ = '-e'^^'r^,, = -/^^%kir^^ = '^{f\-m . (13) 

The antisymmetric tensor eijk in Eq. (fT^ is contracted with spatial components of T^" 
and V'^'^ in Eq. Q. The corresponding quantities T^^eijk and V^^eijk transform, respectively, 
as an axial and a polar vector under the extended rotational group. Therefore, in the non- 
relativistic limit T^^eijk can be absorbed into A^^k and V^^eijk into V^^k- Similarly, taking 
the Levi-Civita structure e^'^^^ejki from Eq. (fT^ we conclude that T-^flie^^-'^^jki = —2T^fii 
transforms as a polar vector and V^ Qie^^^^ejki = as an axial vector. They can be 

absorbed into, respectively, V^^i and A^^i. With the above considerations in mind Eq. ^ 
can be effectively rewritten in the non-relativistic limit as 

= {S^ + GsNN + (v; + GaNi^i^N , (14) 

where S^^, V^, V^, and A^^ have been redefined to absorb the contributions from the vector 
and tensor interactions. Moreover, since the general arguments used for Vn and also 
apply to and (with the suppression |g|/M^ ~ 10^^ in this case) we can safely neglect 
the latter. 

The scalar density S^{x) in Eq. p4|) can be written as x'^X- For a WIMP of spin J^, x is a 
(2J^ + l)-component non-relativistic spinor {e.g. in the MSSM, x is a two-component Weyl 
spinor). The spatial components of the axial- vector density A''^{x) in the non-relativistic 
limit are proportional to the spin density x^'S'^X- Finally, we arrive at the following form of 
the WIMP-nucleon interaction Lagrangian in the non-relativistic limit, 

= 4:fNxhvNVN + lQV2GFaNX^S^XVNSNVN + O [j^-^ ' (1^) 
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where t^tv is the two-component Weyl spinor for the nucleon (initial and final state spinors 
may be different), N = n,p, S^- = c?/2, and the notation for the couplings has been adjusted 
to match that of Eq. (P). In this form, Eq. (|TH|l is valid for non-relativistic point-like WIMPs 
of arbitrary spin. 

The operators neglected in Eq. (fTH|l are suppressed by a factor of 0.01\/A relative to the 
leading order terms [see Eq. ©)]• Strictly speaking, in neglecting such terms we made an 
imphcit assumption that the couplings /tv and Gfcln in Eq- ffT^ are not suppressed relative 
to the other couplings in Eq. such as Gp. If it turned out that Gp > IOOx/A/at then one 
would have to consider the N'y^N operator even though it is formally suppressed in the non- 
relativistic limit. In this paper we ignore this possibility, and modulo the aforementioned 
assumption our analysis is model- independent. 



A. Some Examples 

Below we consider some examples of theories that have WIMPS with spin-0, spin-1/2 
(Dirac particle), and spin-1 and explicitly show that in each case the WIMP-nucleon inter- 
action reduces to Eq. (fT3j) . 



1. Spin-0 

It is possible that in a supersymmetric theory the lightest SUSY particle is the sneutrino 
z/, the scalar partner of the neutrino. In the non-relativistic limit, such a particle interacts 
with a nucleon as P], 

Cj) = —\/2GfJ^Jn,^i , 

= (t^ - 2Qn sin^ Ow) NYN , (16) 

where Pi^,p'c^ are the initial and final momentum of u, and Qn are the isospin and the 
electric charge of the nucleon, and 6^ is the Weinberg angle. For a slow-moving Dirac 
fermion of mass M the leading terms in the v/c expansion of the u- and f -spinors are (in 
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the Weyl representation), 




where p is the fermion momentum and rj is a Weyl spinor. Using the above expressions one 
can exphcitly verify that only the time component of the nucleon vector current survives in 
the limit q ^ 0. Therefore, one is left with a purely SI interaction, which is a special case 
of Eq. (USD with /jv = -V2Gf {Tif - 2Qn sin^ 6^ and oat = 0. 



2. Spin- 1/2 

In theories with extra dimensions Kaluza-Klein (KK) excitations of the SM particles can 
produce viable dark matter candidates. As an example, consider a Dirac KK neutrino, which 
interacts with nucleons via Z° exchange as 

= -V2Gf - 2Q^ sin^ Ow) ^kk1^.^kkNi^N . (18) 

Just like in the previous example, only time components of currents remain in the non- 
relativistic limit, and we again obtain Eq. ^ with Jn = -\FiGf {t% - 2Qns\\^ Q^) and 

Oat = 0. 



3. Spin-1 

Kaluza-Klein excitations of the SM gauge bosons could constitute cold dark matter [3| ■ 
If, for example, dark matter is composed of excitations of the hypercharge gauge boson Bi 
the WIMP-quark interaction has the form: 



7g 



- + 7,) Bl^B.^qq - ta,BlB,,e''f^^Pq^^^,q , 



2 cos2 e 
E. 



2 

m 1 
1l 



m 



Bi 



^2cos2^ 

^2 



m 



4 cos2 6 m\ ' 



(19) 



(20) 
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where Yq is the hypercharge of the quark with mass ruq and energy Eg, rriqi is the mass of 
the quark's first KK excitation, and rrih is the Higgs boson mass. The first term in Eq. (|TI?]) 
corresponds to the SI interaction in Eq. (fT3|) : 



niq ^ rUq 
= ^ E (/^. + 7.)^/^, (21) 

u,d,s,... 9 



where for various quark flavors can be found e.g. in Refs. [I|,ll3[, and /at is normahzed 
to reproduce Eq. (jl5p . The second term produces a SD WIMP-nucleon interaction: 

ze''''''BlB,,{N\q^kl5q\N) = -4A^ {B,\Sb\B,)NSnN , 

= E «.Af , (22) 

where _ = ie^^'^^ is the matrix of the spin operator for Bi, is a dimensionless 
constant of order unity, and is normahzed to reproduce Eq. A recent analysis gives 

AP = = 0.78 ± 0.02, = = -0.48 ± 0.02, and = A^ = -0.15 ± 0.02 L 



III. WIMP-NUCLEUS CROSS SECTION 

The SI and SD WIMP-proton cross sections at low momentum transfers are easily cal- 
culated from Eq. (fT3j)^. 



SD 128J;^(Jx + 1) Jp( Jp + l)G|'a^ 



MrediM) = , (23) 

where and are the WIMP mass and spin. In our convention, the non-relativistic 
spinors ip for both the WIMP and and the nucleon are normalized as ip\pipx,p ~ ^x,p- '^^^ 
SI and SD WIMP- nucleus cross sections at asymptotically small energies, ctq^^ and ctq^, 
can be expressed in terms of the corresponding WIMP-proton cross section, 

,1 _ Ml^{Mnu 



,2 . , . X r . n 2 



Z+{A-Z)^ 
jp 



SI 
XP 



^ For a non-self-charge-conjugate WIMP, such as heavy Dirac neutrino, the WIMP-proton cross section 
must be multipUed by an extra factor of 1/4. This modification does not affect Eq. H24|) . 
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0", 



SD 
OxN 



Ml,{M^^,)4{J + l) 



3J 



{Sp) + {S„)- 



cr. 



SD 

XP 



(24) 



Here, A and Z are the mass number and the charge of the nucleus with spin J, and fp{fn) 
and ap(a„) are the SI and SD WIMP-proton(neutron) couphngs, respectively. The quantities 
(Sp) and (Sn) are the average spins of the proton and the neutron in the nucleus. At finite 
momentum transfer one must average the one-nucleon operators from Eq. over the given 
nucleus using some nuclear structure model. Including the WIMP velocity distribution (see 
Ref. Q]), one then obtains for the elastic scattering rate of WIMPs on the nucleus (per unit 
detector mass), 



Px 



dR 



dE 4veM^MI,{M^ 



erf 



+ Ve 



— erf 



+ Ve 



SI p2 (Tr\j_^SD Sa{E) 



Sa{0) 



EMr, 



(25) 



where is the local Galactic halo density ^, is the recoil energy of the nucleus with mass 
Mnuci and ve and vo are, respectively, velocities of the Earth and the Sun in the Galactic 
frame. 

The form-factors Fsi{E) and Sa{E) ^^j^' depend on the nuclear structure. The SI 

3ji(gi?i^ 



form-factor can be well approximated by 

Fsi{E) 



qRi 



Ri 



'^MnucE 



1.44^2/3 _ 5^2 

1 fm , 



(26) 



where ji{x) is the spherical Bessel function. Unfortunately, there is no such universal ex- 
pression for Sa{E). It can be parameterized in terms of three nucleus-dependent functions: 

1 



Sa{E) 



{ap + an) Soo{E) + (op - a„) Sn{E) + (a - a„)5'oi(^) 



(27) 



where Sij{E) for most nuclei used in WIMP searches can be found in Ref. 17 1. 



We take Px=0.3 GeV/cm^ 1^. 
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It is apparent from Eqs. ()24|) and ()25|) that the detection rate for non-relativistic WIMPs 
depends on only five parameters (in addition to nucleus- dependent constants). We choose 
these parameters to be a^p, cr^p, fn/ fp, (^n/f^p, and M^. Various extensions of the SM 

Snerally reduce the number of parameters {e.g., in the MSSM one normally has fn/ fp ~ 1 
). In our analysis we do not impose such restrictions. 

IV. EXPERIMENTAL DATA 
A. DAMA Nal 

In utilizing the datapublished by the DAMA collaboration in Ref. we adopt the same 
approach as in Refs. |ll|, Q- particular, we define the quantity, 

.2 / ,, «,,„\2- 



( Qih _ riei'pX / qth _ qexp\ ' 



(28) 



where S^^f^^^'^ and S^l^'}'^^^^ are, respectively, the theoretical (experimental) average value and 
the annual-modulation amplitude of the detection rate in the ith energy bin. Theoretical 
predictions for these quantities as functions of M^, a^^, cr^p, fn/ fp, and an/a^ can be 
obtained using Eq. (|^. In this work, we used the experimental values for S^^^^^^- given in 
the first two columns of Table 1 in Ref. 0]. The DAMA preferred region, shown in Fig. 4a) 



of Ref. p, is well reproduced for k ^100. 



B. DAMA i29xe 

The latest results on WIMP searches with a liquid xenon target by the DAMA collabo- 
ration have appeared in Ref. We take the limits on counts per detector per unit mass 
)er day (detector rate unit, dru) for various energy bins from Fig. 4 and Table I of Ref. 



IG I . In calculating the limits on the WIMP cross sections we follow the approach of Ref. 
19|. We take the central values for the dru''s in all energy bins to be zero and the 90% 
CL upper limits on the rfr-u's to be equal to 1.3 times the total error bar on the dru's. We 

n 

verified that the 90% CL upper limits appearing in the last column of Table I in Ref. [10| 
are reproduced in this way with the error bars from the second column of the same table. 
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In this approach, the 90% CL upper hmits on the dru in each energy bin result in an 
upper hmit on the WIMP-nucleon cross section a^^^{k) with the help of Eg. () 25 |) . The 
combined upper limit from all energy bins is obtained using Eq. (15) of Ref. jl9l |. 

1 V . ^7—. . (29) 



We verified that our calculation reproduces the exclusion curves shown in Figs. 6(a,b) of 



I — I 



Ret. m for both the SI and the SD cases. 



C. EDELWEISS 

We used Ref. |3] to incorporate the latest results from the EDELWEISS experiment that 
used a heat-and-ionization cryogenic Ge detector. In order to extract the 90% CL limit on 
the average expected number of events N from the null experimental result we used the 
Bayesian approach described in Ref. [2^. For our analysis, we adopt the following prior 
distribution for A^, 

, , f for AT < , 
^W={l foriV>0; 

which leads (after normalization) to a PDF for of the form po{N) = . Therefore, the 
probability for to be below some value A"o is: 

/•No 

PoiN <No)= PoiN)dN = 1 - e-^° . (31) 
Jo 

Equating the above probability to the confidence level of the EDELWEISS constraint (90%), 
we obtain A"o = In 10 ~ 2.3. 

The predicted number of events in the EDELWEISS detector for fixed values of WIMP 
parameters can be related to the detection rate by integrating Eq. (f^3j) over the range of 
energies accepted by EDELWEISS (20 keV to 64 keV) and multiplying by the total effective 
exposure of 11.7 kg-days Qj. In order to properly normalize the detection rate we multiply it 
by a factor Ce- This factor is adjusted such that the following predictions given in Ref. [8| for 
numbers of events are reproduced for the pure SI case: for M^=44 GeV and a^p = 5.4 x 10~^ 
pb one has A^ = 6.2, and for M^=52 GeV and = 7.2 x 10"^ pb one has A^ = 9.8. Note 
that since there is only one coefficient to determine, the second data point is redundant, and 
can thus be used for a consistency check. We find that for Ce ~0.76 both test points and 
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the EDELWEISS exclusion curve given in Fig. 5 of Ref. P] are well reproduced. In all the 
Figures shown below, the constraints from EDELWEISS data are obtained by setting the 
predicted number of events in the detector equal to A^o = In 10. 



D. ZEPLIN-1 

The ZEPLIN-1 experiment (see Ref. ^j) used liquid Xenon as a medium and relied 
on pulse-shape discrimination analysis. Unfortunately, at present there is no publication 
available containing detailed experimental results. The absence of such data prevented us 
from treating ZEPLIN-1 results in the same manner as those from the DAMA experiment 
with liquid ^^^Xe (see above). Instead, we chose an approach similar to the one we used 
for the EDELWEISS experiment. We deduced the parameters necessary for the calculation 
from Ref. j^. In particular, we chose the visible energy to be between 4 keV and 30 keV, 
with quenching factor gxe(Zeplin) = 0.2. These parameters correspond to recoil energies 
between 20 keV and 150 keV. In order to calibrate our calculation we introduced an effective 
statistics factor, Cz- We found that the SI 90% CL ZEPLIN-1 limit shown in Ref. Q is 
exceptionally well reproduced for ^ 8 kg-days. 

Due to the absence of published results, the ZEPLIN-1 limit is the most uncertain input 
in our calculation. However, the results of our analysis are robust and are unlikely to change 
when the full data is properly included. Indeed, no change in position of the ZEPLIN-1 
curve in Fig. |21will move the upper limit on the WIMP mass above about 25 GeV, which 
is still in conflict with the present lower limit on the MSSM lightest-neutralino mass of 37 
GeV. Similarly, even if the ZEPLIN-1 constraint is entirely removed from the analysis one 
would still require fn/ fp < to achieve even marginal agreement among all data for the 
predominantly SI case (see explanation to Fig. As shown below fn/ fp < would be 
very unusual in the MSSM. 



E. Energetic neutrino searches 

A promising method for indirect detection of WIMPs is the search for neutrinos from 
WIMP annihilation in the Sun and/or the Earth. Such neutrinos produce upward muons 
in the Earth via charged-current interaction. Measurement of (or constraint on) the flux of 
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such muons indirectly constrains the annihilation rate of WIMPs, which is related to the 
WIMP capture rate by the Sun and/or the Earth. A detailed review of such indirect WIMP 
detection with extensive list of references on the subject can be found in Ref. Il|. In this 
work, we use the limits on the fluxes of neutrino-induced upward muons near the surface of 

n 

the Earth obtained by Super-Kamiokande 

r^(sun) ~ r^(Earth) < 10"^ m"^ year"^ . (32) 

Limits from other experiments, such as 1MB, Baksan, MACRO, AMANDA, etc., are very 
similar (see Refs. |l[ 1ll| and references therein). The upward muon flux is related to the 
WIMP-proton cross section as 



r- = 1.96xlO-dta„h^^-j«M,)f(M,)^^j ^^^^ 
r^°(sun) = 1.6 X 10-nanh= (1) {(Af,)S(A/,/M,) {■^) (^ 



SD 

> 1 2 



(33) 



where d = 3.3 x 10* m~^ year~^ for the Sun and d = 1.7 x 10* m^^ year^^ for the Earth, and 
T is the time scale for equilibration between WIMP capture and WIMP annihilation. For 
the Sun, we can take tanh(t/r) ^ 1, whereas for the Earth it is likely substantially smaller 
than unity Q. The function ^(M) is given by Eq. (9.54), and S{x) by Eq. (9.21) in Ref. Q. 
The function f'{M) is the generalization of the function f{M) given by Eq. (9.28) of the 
same reference, 



f (M) =Y.lAS{MjMi^jF,{M^ 



nuc X 



_|_ A _ f_u_ 



(34) 



M2(Mx + M„ 

where the sum runs over the nuclei with mass M^^^,, charge Zi, and atomic number Ai. The 
quantities and (pi are given in Tables 8 and 9 of Ref. flj. The form-factor suppression 
Fi{M^) is only important for iron in the Earth, where it is given by Eq. (9.23) of Ref. Q]. 

For ^ Mp the muon flux T^^ is roughly independent of M^, which leads to the 
constraint a'^^ < 1.2 x 10"^/.^ pb. In the MSSM one generally has C, >0.03 1], and we 
obtain a^^ < 4 x 10^'^ pb. The constraints on a^p from energetic solar neutrinos and 
on (T^^p from terrestrial neutrinos are too weak to be interesting. On the other hand, the 
constraint on a^p from energetic neutrinos originating in the Earth is potentially non-trivial. 
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For the case fn/ fp ~ —0.76 and ~ 50 GeV (see Section IV B|) one obtains 



^XP < .7^2^/ ^ Pb • (35) 



3.6 X 10 



~ etanh'[t/r, 



As shown below, one needs a^i ~ 0.0035 pb for this case. In addition, one can deduce from 



XP 

Ref. |2l|^ 

f-V = 2-4 X 10V'(50 GeV) J-^^^^ ^ , (36) 

\r©y 10 '^^ cm'^s ^ Ipb 

where is the WIMP annihilation cross section in the limit of zero relative WIMP ve- 
locity V. In order to have a WIMP relic density of order unity we must have a^v ~ 
10~26 cm^ s-i and for fn/ fp = -0.76 one has /'(50 GeV) = 0.85. Therefore, 

{t/T(sY ^ 2 X 10^ (a^^/lpb) ^ 8, or tanh^(t/re) ^ 1. Now, Eq. §^ can be simply 
rewritten as a constraint on ^ < 0.001 for fn/ fp = —0.76. For the largest value ^ = 0.001 
allowed in this scenario the energetic neutrino constraint on the WIMP-proton SD cross 
section is a^^ < 0.12 pb. 

Can the constraints coming from non-observation of upward muons from energetic so- 
lar neutrinos be evaded? Below we consider two possibilities. In the first case, the 
MSSM WIMPs predominantly decay into light fermions. Because the annihilation rate 
T{xx — ^ f^f ) is proportional to the fermion mass squared m'j, direct annihilation into 
neutrinos is virtually impossible, and energetic neutrinos appear in the decay chain of the 
initial annihilation products. The total flux of the neutrinos originating from the branch 
T{xx ~^ f^f ) inherits the suppression by m'j, and in the case where only light fermions 
appear during annihilation the flux may be orders of magnitude below the conventional esti- 
mates Effectively, in this case ^(M^) is substantially smaller than estimated in Ref. 
which weakens constraints on cr^p^^- In the second case, the WIMP is not identical to the 
anti-WIMP, and only WIMP-anti-WIMP annihilation is allowed. Although direct annihi- 
lation into neutrinos may be possible (leading to stronger energetic neutrino signals), the 
annihilation rate may still be significantly suppressed in the presence of significant WIMP- 
anti-WIMP asymmetry. 



1. WIMPs only decay into light fermions 



As pointed out in Ref. ll| the energetic-neutrino constraint could be evaded if the 
WIMPs annihilated to uu, dd, ss, e+e~, and/or pairs but not cc, bb, nor rf pairs. 
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Such a situation can in principle be achieved in the MSSM by fine tuning the sfermion 
masses and the WIMP composition. Let us write the hghtest neutrahno field as 

X = ZiB + + Z3H1 + Z^H2 , (37) 

where Zi are constants subject to J2 \Zi\^ = 1, B and W3 are superpartners of the B and W3 

gauge bosons, and Hi^2 are superpartners of the neutral Higgs bosons. Examination of the 

general expressions for the axial-vector and scalar neutralino-fermion couplings (Eq. (3.6) 

in Ref. 13|) shows that if IZ3I = IZ4I and Z2 = tan -2^1 the scalar coupling vanishes and 

the axial coupling a/ is inversely proportional to the fermion superpartner mass squared, 

1 . , 

«/ ~ T72 (38) 

/ 

where Mj is the mass of the superpartner of the fermion /. Choosing the superpartner mass 
to be very large for the charm and bottom quarks, as well as for the r lepton, one can force 
the neutralinos to annihilate into the light quarks and leptons only. However, choosing such 
fiavor non-universal masses for the scalars may be problematic in the MSSM. Indeed, one 
must obey the existin g st ringent constraints on the size of fiavor- changing neutral currents 
(FCNC)(see e.g. Ref. for a review of FCNC constraints on the MSSM spectrum). Since 
evading energetic- neutrino constraints in a way we just described requires significant fiavor 
non-universality among certain entries of the squark mass matrices, whereas smallness of 
the FCNC prefers the opposite, one may expect that fine tuning would be required to make 
any such scenario phenomenologically viable. We note that even for non-MSSM WIMPs 
fiavor non-universality of the WIMP couplings is needed to evade the energetic-neutrino 
constraint in the described manner, and experimental limits on FCNC are likely to present 
complications for any WIMP candidate. 

2. The WIMP is not identical to anti- WIMP 

Another way to evade the energetic-neutrino constraint is to consider WIMPs that are 
not identical to their antiparticles. Indeed, for the total number density of the dark matter 
particles Nq one has 

no + no = No , 

no -no = W , (39) 
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where Uq and no are the WIMP and anti-WIMP densities in the vicinity of the Solar system, 
and W is the "WIMP number" density, which may be generated in the presence of some 
T-, C- and CP-violating physics. The (anti-)WIMP densities n and n inside the solar core 
obey the equations 

n = AriQ — Bnn — T^sc , 

n = Auq — Bnn — Tesc , (40) 

where the first and the second terms on the RHS of both equations represent, respectively, 
the capture rate of WIMPs incident on the Sun and the annihilation of WIMPs and anti- 
WIMPs in the solar core. The constants A and A are proportional, respectively, to the 
WIMP-proton and anti-WIMP-proton scattering cross sections. The last terms represent 
the losses of WIMPs due to evaporation (see e.g. Ref. 

Consider the case where uq ^ tiq. In the stationary state n = 0, and we have for the 
annihilation rate Tann = Ano. We make a reasonable assumption that annihilation is the 
dominant mechanism for anti-WIMP loss if -C because the evaporation rate is suppressed 
by a small Boltzmann-like factor ^-^e/kTw _ Here, is the escape energy and Tw is the 
effective WIMP temperature. On the other hand, for neutralinos in the MSSM we have 
{Tesc can be neglected here as well). 

In the stationary state this would give for the annihilation rate Pq^!^"^*^ = A^^^^'^^ Nq. We 
find that 

^ A flQ MSSM ^ ^ -pMSSM f 1 _ ^\ (Ac^\ 

''-ann A MSSM ann at ) \^^) 



ann 



j^MSSM ]\J-^ ann 2 A^'^ ^ ^ \ Nq / 

We generally expect A ~ A, and since the WIMP-proton cross section is constrained by 
DAMA we expect A ~ A*^^"^^^ Therefore, we conclude that T^nn ~ r*^^^*^(l - W/Nq). 
The flux of upward muons in the Super-Kamiokande detector is proportional to the WIMP 
annihilation rate: P^ ~ C,{M^)rann- Therefore, we obtain, 

e(Mj _ W\ ^MssM _ (43) 



^MssM^M^^ y pjj 

Since the WIMP is not identical to its anti-WIMP, direct WIMP annihilation into neutrinos 
may be possible, which may yield E,{M-^ / E,^^^^^'^ [M^) > 1. On the other hand, for 1 — 
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W/Nq <^ 1, the muon detection rate in the Sun can still be significantly below the one 
predicted by Eq. (5) in Ref. 11[, and the energetic- neutrino constraint may not apply. In 



this work, we first perform the analysis without this constraint, and then add it later on. 



V. ANALYSIS 



A. The energetic-neutrino constraints are not applicable 

Let us first consider the constraints on the WIMP parameter space without the energetic- 
neutrino bounds. In this case an agreement between all direct search experiments can be 
achieved for a predominantly SD WIMP-nucleon interaction for a wide range of WIMP 
masses aud coupl.ugB. A similar couclus.ou was or,g„>ally obtained ■„ Ref. U In F,g. □ 
we plot the allowed region in a^p vs. plane for fp/ fn = 1, CLn/dp = and =50 GeV. 
The shaded region is allowed by all direct search experiments described in Section IIVI Since 
a substantial SD interaction is required, the WIMP cannot be a scalar particle in this case. 



B. The energetic-neutrino constraints are applicable 

We now include the energetic- neutrino constraint. As is well known, for the usual case 
fn/ /p ~ 1 the results from all direct search experiments cannot be reconciled for the pure SI 
case at about 3a level, and some SD interaction appears necessary. With the WIMP-proton 
SD cross section limited from above by the energetic-neutrino searches, a relatively large 
WIMP-neutron SD cross section is needed to accommodate the DAMA/Nal result. The situ- 
ation is similar to the one described in Ref. jll| for the pure WIMP-neutron interaction, and 
the allowed region for this case is shown in Fig. |21 For this figure, we set the WIMP-proton 
SD cross section to zero thereby automatically satisfying the energetic- neutrino constraint. 
In this scenario, there is an upper limit on the WIMP mass. The strongest limit, about 18 
GeV, will likely be provided by the ZEPLIN-1 experiment. We found that this limit does 
not significantly change even if some SI interaction is allowed. If the ZEPLIN-1 result is not 
included the upper limit on the WIMP mass increases to about 25 GeV. 

It is interesting to note that MSSM neutralinos with mass below 37 GeV are excluded by 
direct colhder searches . Although the analysis of such searches is not completely general 
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FIG. 1: Constraints in the a^p vs. cr^^ plane due to various direct WIMP searches provided that 



XP 



the energetic-neutrino constraint is not included. In this plot, fn/fp = 1, an/ dp = 0, and 



50 



GeV. The shaded area is the 3a allowed region consistent with the annual modulation observed by 
the DAMA collaboration in their Nal detector j^]. The region below the dashed curve is the region 
allowed by the DAMA ^^^Xe experiment 



jlOIJ . The region below the dotted line is that allowed by 
the EDELWEISS experiment and region below the dash-dotted line is allowed by ZEPLIN-I 
experiment 



since it assumes gaugino and sfermion mass unification at the GUT scale j20[, evading this 
limit, if at all possible, would require fine tuning of the MSSM parameters. 

Up to this point, we always maintained the condition fn/fp = 1- It turns out that 
relaxing this constraint allows one to achieve marginal agreement among all data for in 
the vicinity of 50 GeV. However, with the ZEPLIN-1 result included this occurs for only a 
very narrow range of fn/fp- The possibility arises because the neutron-to-proton ratio differs 
slightly from one nucleus to the other. Specifically, if —0.77 < fn/fp ^ —0.75, then the 
WIMP coupling to the neutron in Xe cancels that from the proton, allowing a null result 
in ZEPLIN-1 to be consistent with a modulation in DAMA. If ZEPLIN-1 is not considered 
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FIG. 2: Constraints in the cr^„ (pb) vs. (GeV) plane due to various direct WIMP searches. To 
ihustrate, we took ap/an=0 (pure WIMP-neutron interaction). This choice automatically ensures 
that the energetic-neutrino constraints are satisfied. Notation is the same as in Fig. ^ 



the allowed range is increased but the consistency range for /„/ fp is still very narrow. The 
situation for fn/fp = —0.76 is shown in Fig. |21 The small region allowed by all the data 
is centered around a^p ~0.0035 pb. Note that, due to significant cancellation between the 
WIMP-neutron and the WIMP-proton scattering amplitudes, the individual WIMP-nucleon 
cross sections are about 10^ times larger than in the case where fn/fp = 1- The inverse of 
this number, 0.01, refiects the amount of fine tuning required for this solution to work. 

From the standpoint of the MSSM, this solution is in fact "doubly fine-tuned" . In addition 
to ,„aki„g sure that /„//, . carefully selected to fit all the data oue has to tuue the model of 
SUSY breaking to even obtain fn/fp < (normally, one has /„ ~ fp in the MSSM [Ij). It is 
interesting to ask, therefore, whether a significant deviation from the approximate equality 
fn/ /p ~ 1 is at all possible within the MSSM. In general, we have. 



fN = fo + fih , 



(44) 
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FIG. 3: Constraints in the (T^p(pb) vs. (7^„ (pb) plane due to various direct WIMP searches. We 
took /n//p=-0.76. Notation is the same as in Fig. ^ The energetic-neutrino constraint for this 
case, < 0.12 pb (see Section II V E() . is satisfied. 



where N = p,n and is the usual Pauli matrix in the strong isospin space. The Feynman 
diagrams and detailed expressions for /p „ can be found in Ref. |]|. Since only the up and 
down quarks have non-zero strong isospin, the isovector part of Jn is entirely due to coupling 
of the neutralino bilinear to isovector operators built from the up and down quarks. The 
isoscalar part arises due to coupling to various isoscalar operators built from the up and 
down quarks, the strange and all heavy quarks, and the gluons. The scalar coupling of the 
neutralino to a quark has the form 

where Gp is the Fermi constant, Ag and Bg are dimensionless constants, rrig is the quark 
mass, Mg is the mass of the superpartner of q, Mh is the lightest Higgs-boson mass, and 



^ We assume that the squark mass matrices are diagonal in the flavor space. 
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is the neutralino mass. We can now write the couphngs /p^„ in the form 

fp,n 



a, 



fnp,n [Itu fu + frd fd + It's fs + 



(46) 



where "■ ■ ■" stand for both the isovector terms of order 1/M| and all remaining isoscalar 
contributions. Here, mi^if^^ = {N\mgqq\N). Ignoring the remaining terms we obtain (see 



Table 6 in Ref. 



fl 


< 


fuifxu fru) + fdifrd frd) 




-0.014/, + 0.025/rf 


fo 




2fTsfs + WL + f^u) + fdifrd + frd) 




/, + 0.15A + 0.27/d 



(47) 

) JdyjTd^ JTd) js -r v.Loju jd 
In most models of SUSY breaking /„ ~ fg, and although the above range is nothing more 
than an estimate it at least shows that |/i//o| > 1 (which would lead to a negative ratio 
fn/ fp) is generically disfavored in the MSSM. However, it is possible in principle to fine-tune 
the parameters to obtain |/i//o| > 1. One scenario arises when |/s| ^ \fu,d\- This can be 
achieved if Mg is very large for all squarks except for u and d, all physical Higgs bosons, 
except the lightest one, are very heavy, and coupling between the WIMP and the lightest 
Higgs boson is tuned to zero. In this case, one finds (we use the MSSM Feynman rules given 
in Ref. [2^ and bring the notation in correspondence with Ref. jl^), 

,2 



fu 

fd 
fs 



9 



-Z4 {Z2 — tan 6^14/ Zi) 



AMwMi sin p 
2 

Z3 {Z2 - tanOwZi) 



9^ 



AMwMj cos (3 
. 



(48) 



Here, g is the SU{2)l gauge coupling, and tan/5 is the ratio of the expectation values for 
the two Higgs boson doublets in the MSSM. We find 



fo 



0.014 + 0.25 tan^/? 



2 /,Af| 

d 



-0.15 + 0.27tan2/5 



(49) 



If M^/Mu ~ 1.3 tan /3 then one can have |/i//o| > 1, which may lead to fn/ fp < 0. At 
present, tan/? > 3 is favored by precision , one needs mild hierarchy 

between the up- and down-squark masses to obtain fn/ fp < 0: > 3.9Mu. 

There might be other scenarios yielding /„ substantially different from fp. However, they 
will all have to share the same property: the first generation of squarks must be singled out 
from the rest to enhance the up- and down-quark contributions to /o,i. Reconciling such 
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significant flavor non-universality of tlie squark flavors witli stringent constraints from FCNC 
will generally require flne-tuning. In this sense, having signiflcantly different values for 
and fp is unnatural, although potentially possible, in the MSSM. There is an additional 
complication, however. For this scenario to work, the WIMP-proton scattering cross section 
must be roughly 100 times larger than for the case fn ~ fp- On the other hand, one 
must have ^ < 0.001 in this case (requiring Z2 — tanOwZi <^ 1, see Section |IVE1|) . which 
suppresses fu,d according to Eq. (PH|) . Partly due to this suppression, we have found no 
scenario where sufficient enhancement of fp^n would occur in the MSSM for ^ 50 GeV 
together with fn/ fp ~ —0.76. 

We point out that fn/ fp ~ —0.76 does not hold for any of the dark-matter candidates we 
considered. For scalar neutrinos and KK neutrinos one finds (see Subsection III A|l /„/ fp = 
T^/{Tp — 2QpSm^ 9iy) ~ —10. A detailed study of this ratio for KK excitation of the 
hypercharge gauge boson shows that it is at most few percent away from unity. In summary, 
although there exists a small phenomenologically allowed region of parameter space where 
predominantly SI WIMP-nucleon interactions can marginally account for all data available 
on WIMP searches this region appears to be out of reach for all of the WIMP candidates 
we considered. 

VI. CONCLUSIONS 

In this work we performed a generalized analysis of the dark-matter detection experiments 
listed in Section HVl Our analysis is formulated in terms of the WIMP mass M^, the SI and 
SD WIMP-proton cross sections a^p'^^, and coupling ratios fn/ fp and an/ap. We found 
several regions in this parameter space that allow for agreement among all data. 

• If the energetic-neutrino constraint does not apply it is possible to reconcile all data 
with a predominantly SD WIMP-proton interaction. This scenario can in principle 
be tested in a direct dark matter detection experiment with an odd-proton target, 
such as e.g. ^^F (26]. Evading the energetic-neutrino constraint is problematic in 
the MSSM because it generally requires highly fiavor non-universal squark masses. 
Another possibility is a non-self-charge-conjugate WIMP, which would indicate physics 
other than the MSSM. In order for this possibility to be realistic, such physics must 
accommodate WIMP-number-violating operators and possess a sufficient amount of 
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C and CP violation to generate a fractional WIMP asymmetry near unity. 

• If the energetic-neutrino constraint applies then the SD WIMP-proton cross section is 
constrained to be a^p < 1.2 x 10~^/^ pb, and for a wide range of fn/ fp, a SD WIMP- 
neutron cross section cr^^ > 30 pb is required. In addition, is constrained to 
be below 18 GeV, which is inconsistent with the present lower limit on the lightest 
neutralino mass within most models of SUSY breaking [2^ . 

• A heavier WIMP (about 50 GeV) is allowed for a narrow range —0.77 < fn/ fp ^ 
—0.75 with essentially no SD interaction present. In addition, one must have 
a^p ~ 0.0035 pb. None of the WIMP candidates we considered could satisfy both 
requirements at the same time. 

We see that although it is possible to reconcile the direct and indirect dark-matter search 
experiments listed in Section IIVI the resulting parameters are in general unnatural for the 
MSSM with most models of SUSY breaking that appear in the literature. Either significant 
flavor non-universality among the SUSY breaking parameters {e.g. to evade the energetic- 
neutrino constraint) or some delicate relationships among them (to avoid the lower limit on 
the lightest neutralino mass) appear necessary. Generating the WIMP parameters in any 
of the allowed regions may be problematic in a theory with non-MSSM WIMPs as well. 
For instance, evading the energetic-neutrino constraint in a theory with non-self-charge- 
conjugate WIMPs may require a significant amount of CP-violation to generate a fractional 
WIMP number close to one. Such CP-violating interactions may be strongly constrained 
by existing limits on e.g. neutron and atomic electric dipole moments [2^. If the energetic- 
neutrino constraint is adopted, one would have to explain how very light WIMPs (M^ < 18 
GeV) with significant WIMP-quark couplings have so far escaped detection. 

On the other hand, if the DAMA result is removed from the analysis the MSSM neutralino 
(as well as other WIMP candidates) will be compatible with the remaining experiments and 
still be capable of being produced with sufficient abundance to account for all dark matter 
I20I. In view of this situation, one can hope that forthcoming direct dark-matter search 
experiments, potentially capable of improving the present sensitivity by several orders of 
magnitude, will be able to conclusively confirm or exclude the results published by DAMA 
based on their observation of annual modulation of the detection rate. 
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